Competition of isoscalar and isovector proton-neutron pairing in nuclei 
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We use shell model techniques in the complete pf shell to study pair correlations in nuclei. 
Particular attention is paid to the competition of isoscalar and isovector proton-neutron pairing 
modes which is investigated in the odd-odd N = Z nucleus 46 V and in the chain of even Fe-isotopes. 
We confirm the dominance of isovector pairing in the ground states. An inspection of the level 
density and pair correlation strength in 46 V, however, shows the increasing relative importance 
of isoscalar correlations with increasing excitation energy. In the Fe-isotopes we find the expected 
strong dependence of the isovector pairing strength on the neutron excess, while the dominant J = 1 
isoscalar pair correlations scale much more gently with neutron number. We demonstrate that the 
isoscalar pair correlations depend strongly on the spin-orbit splitting. 
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I. INTRODUCTION AND MOTIVATION 

Stimulated by progress in many-body theories and the 
advent of novel experimental tools like the radioactive 
ion-beam facilities, pairing in nuclei, and in particular 
proton-neutron (pn ) pairing, has become quite fashion- 
able recently pi rLOfl. In contrast to like-nucleon pairing, 
pn pairing may exist in two distinct varieties - isovector 
and isoscalar pairing. Isovector pn pairing has been found 
to decrease fast with increasing neutron excess (or more 
precisely with increasing \N — Z\, where N, Z are the 
numbers of neutrons and protons, respectively) within an 
isotope chain |[. Thus, it is obvious that pn pairing is 
relatively more important in N = Z nuclei. While even- 
even N = Z nuclei have ground state isospin T = and 
isospin symmetry forces isovector pairing to be identical 
in all three T — 1 pairing channels (proton-proton (pp), 
neutron-neutron (nn), and proton-neutron), this symme- 
try is broken in most odd-odd N = Z nuclei with mass 
number A > 40 which have T = 1 ground states. Indeed, 
it is observed that isovector pn pairing is the dominating 
pairing mode in these nuclei. On the other hand, odd- 
odd N — Z nuclei with A < 40 have ground state isospin 
T = (with the only exception of 34 C1) pointing to a 
competition between the two different pn pairing modes 
and its change with mass number is an obvious source 
of interest. The interplay of isovector and isoscalar pn 
pairing has been studied in various contexts: 

• It has long been known that N = Z nuclei gain 
extra binding compared to the smooth trend de- 
fined by their neighbors. This extra piece, called 
the Wigner energy, is stronger for even-even than 
for odd-odd nuclei. Recent studies of the Wigner 
energy tried to relate its origin to the pn pairing 

• Odd-odd N = Z nuclei are the only ones with two 
nearly degenerate low-lying isospin states (T = 



and T = 1). This fact and the change of the ground 
state isospin from T — to T = 1 at A = 40 has 
been explained by a near balance of isovector pair- 
ing and symmetry energy flf][| . 

• While the ground state isospin of odd-odd N = Z 
nuclei with A > 40 (except for 58 Cu) is T = 1, it is 
experimentally found that the density of T = lev- 
els at modest excitation energies is larger than the 
density of T = 1 states. For example, for 74 Rb a 
transition from isovector to isoscalar dominance has 
recently been observed with increasing rotational 
frequency p"2| . This transition has been traced 
back to a competition between isovector J = pn 
pairing and isoscalar aligned pn pairing (J = 9), 
where the latter becomes obviously more important 
at larger angular momenta Jt3[ . 

The interplay and competition of isovector and 
isoscalar pairing is also the theme of this paper, where we 
study and discuss two different aspects of this topic. At 
first, we perform detailed studies of the spectrum of V. 
Special emphasis is paid to the comparison of the T = 
and T = 1 parts of the spectra and their dependence on 
the pair correlation strength. In particular, we investi- 
gate the temperature dependence of the various pairing 
modes tracing possible differences in the behavior of the 
isovector and isoscalar pn pairing. 

Isovector pairing in the dominating J — channel in- 
volves coupling of states in time-reversed components of 
the same orbitals. Very often isovector pairing is ef- 
fectively given by one single matrix element describing 
the J = pairing of the valence nucleons in one partly 
filled orbital (the fan orbital in the nuclei studied in 
this paper). In contrast, isoscalar pn pairing, even its 
L = component, can involve not only pairing within 
the same orbital, but also between spin-orbit partners. 
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Thus isoscalar pn pairing involves several scales: the cou- 
pling matrix element (s) in the same orbital, the usually 
strong (e.g. J — 1) matrix element between spin-orbit 
partners, and finally the spin-orbit splitting which can 
be understood as a penalty which reduces this otherwise 
very favorable correlation. To shed some light on the 
importance of isovector versus isoscalar pairing, and in 
particular, on the role played by the spin-orbit splitting 
we study pair correlations within the even Fe-isotopes 
50-56p e com p ar i n g the results of "realistic calculations" 
to those in which we have artificially reduced the spin- 
orbit splitting. 

Our studies have become feasible due to recent progress 
made in solving the interacting shell model, which is the 
method of choice to investigate correlations among nu- 
cleons. All calculations are performed with the realistic 
KB3 interaction [[blj in the complete pf shell. Modern 
diagonalization codes allow now studies in model spaces 
whose dimensions have been considered untractable only 
a few years ago. We use such a state-of-the-art diago- 
nalization code p| in our calculation of 46 V and hence 
have the possibility to study the spectrum and its prop- 
erties in a state-by-state approach. About a year ago, 
when we started our investigations of pairing properties 
in the Fe-isotopes, diagonalization codes seemed not to be 
able to handle valence model spaces of these dimensions 
(for very recent progress, the reader is refered to [ff6f ). 
For that reason we have adopted the Shell Model Monte 
Carlo (SMMC) approach for tnese studies. In the 

SMMC the nucleus is described by a canonical ensemble 
at finite temperature. Since the Monte Carlo techniques 
avoid an explicit enumeration of the many-body states, 
they can be used in model spaces far larger than those 
accessible to conventional diagonalization methods. As 
a disadvantage compared to these methods, the SMMC 
approach cannot study detailed spectroscopy and ground 
state properties are achieved as the limit of low tempera- 
tures. To circumvent the "sign problem" encountered in 
the SMMC calculations, we adopt the "g-extrapolation" 
procedure suggested in Ref . |20) . 

Both methods - the diagonalization approach and the 
SMMC - are well documented in the literature and thus 
there is no need for a repetition here. For detailed de- 
scriptions of the methods we refer the reader to . 



II. RESULTS AND DISCUSSION 



First we discuss results obtained for the odd-odd N = 
Z nucleus 46 V by diagonalization in the complete pf 
shell. The calculated spectrum is compared to the ex- 
perimental spectrum [g2| in Fig. [I]. We find overall good 
agreement, although the calculated T = spectrum is at 
slightly too low (by about 400 keV) excitation energy. 

Given the spectrum of the nucleus, we define the in- 
ternal energy at a finite temperature T by (f3 = 1/T) 



U(T) 



l)Eiexp[-{3Ei 



with the partition function 



Z = ^(2J J + l)cxp[-/3£; 4 ] 



(1) 



(2) 



the sums runs over all nuclear levels i labeled by quan- 
tum numbers (Tj, Ji) and excitation energy Ei. Fig. 2 
displays the internal energy and the specific heat C v — 
-£pU(T) as function of temperature T. The striking fea- 
ture is the double-peak structure in the specific heat. 
While the peak around T — 1.5 MeV is trivially related 
to the finite size of our model space (Schottky effect), 
the peak at T = 0.2 MeV is a unique structure of odd- 
odd N = Z nuclei with A > 40 (except for 58 Cu). As 
we pointed out earlier, these nuclei have ground state 
isospins T = 1 and J = 0, and exhibit, like their even- 
even analogs 46 Ti and 46 Cr, a sparse spectrum of T = 1 
levels at low energies reflecting the strong isovector pair- 
ing. However, in odd-odd N — Z nuclei the lowest levels 
with T — and T — 1 are almost degenerate reflecting 
a competition of pairing and symmetry energy. Further- 
more, as can be seen in Fig. 1, odd-odd N = Z nuclei 
exhibit several T = levels at low excitation energy and 
these levels have J > 0. 

The peak in the specific heat at T = 0.2 MeV sig- 
nals the "transition" from the T = 1 ground state to 
the bunch of T = levels (with larger statistical weight) 
at low excitation energy. We note that the peak will be 
somewhat shifted towards higher temperatures if we add 
an attractive component oT 2 , with the isospin operator 
T, to the Hamiltonian to correct for the slight mismatch 
between the T = 1 and T = levels in our calculation. 
In passing we also remark that such a peak is not visible 
in the specific heat for odd-odd N = Z nuclei with T = 
ground states and angular momenta J > 0, where the 
T, J = 1,0 level is at a modest excitation energy and due 
to its small statistical weight does not show up as peak 
in the specific heat. 

To complete our discussion of the 46 V spectrum we 
have plotted the level density in Fig. 3. For clarity we 
have binned the results in 0.5 MeV bins and exhibit them 
for both isospin channels T = and T = 1 separately. As 
is customary, we compare our shell model level densities 
with the backshifted Fermi gas p3[. 



P{U) 



exp 2\ 



12a 1 / 4 C/ 5 / 4 



U = E- A 



(3) 



where A, a are the backshift and level density parame- 
ter, respectively. These two essential parameters of the 
model are determined by global fits to experimental data 
at low energies (around the neutron threshold and be- 
low), leading to the approximate expressions p3[ 



A/8 MeV" 1 

± " - H MeV 

Va a 



(4) 
(5) 
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with the positive (negative) sign for even-even (odd-odd) 
nuclei and A = for odd- A nuclei. We note that the ex- 
perimentally determined value of a is significantly larger 
than the Fermi gas value (a s» A/ 16 MeV -1 ) indicating 
the presence of additional correlations in the low-energy 
levels other than those described by pairing. The calcu- 
lated 46 V level density is reasonably well described in the 
energy regime E > 2.5 MeV by the backshifted Fermi gas 
model with the parameters A = —2.0 MeV and a = 3.5 
MeV -1 . While A agrees well with the empirical value, 
our a is significantly smaller suggesting the presence of 
states outside of our model space in 46 V already at rather 
low excitation energies. The discrepancy in a is less than 
the comparison with Eq. (4) suggests; in Ref. the 
values A = —1.79 MeV and a = 4.9 MeV -1 are recom- 
mended, closer to the ones found here. 

As we have seen already in the discussion of the spec- 
trum, pairing plays an essential role. We will now investi- 
gate this topic in more detail by studying the pair corre- 
lation strength in 46 V as a function of excitation energy. 
We restrict our discussion to s-wave pairing only, hence 
L = 0. As a measure of the pair correlation strength we 
use the operator AfjTt z , defined in LS'-coupling by: 

NjTtz = £ + 1)(2Z' + l)All'] M {ll)A T ^ JM {l'l') 



ll'M 



with the two particle creation operator A^: 

1 r , i i Tt z 



A lsjm( 1 aIb) = 



AD 



a A a B 



LSJM 



(6) 



and a nlms t creates a nucleon of isospin projection t z 
and spin projection s z in the orbital nlm. In our studies 
involving the diagonalization approach the expectation 
value is evaluated state-by-state. In the SMMC approach 
the expectation value refers to a thermal average defined 
as 



(O) 



Tr A exp[-[3H}0 
Z 



(8) 



where Tta is the canonical trace (fixed numbers of neu- 
trons and protons) and H is the many-body Hamiltonian. 

In the limit of large degeneracy (Af) represents the 
number of nucleon pairs with the angular momentum J, 
isospin T and its projection t z . In pf shell nuclei the 
degeneracies are not too large, the operators A, A' are 
not really bosons, and hence (A/j,) can deviate from the 
true pair number expectation value. 

Fig. [| shows the isovector (T — 1, S = 0) and isoscalar 
(T = 0, S = 1) pair correlation strength as function of ex- 
citation energy (averaged in 0.5 MeV wide bins) and for 
T = and T = 1 states in 46 V separately. Isospin sym- 
metry requires (fifj=o,T=i,tz) to be identical for all three 
isovector channels (pp, nn, pn) in T = states, while 
isovector pn pairing can be different from like-particle 



pairing in the T = 1 states (pp and nn pairing strengths 
are still identical). Several observations with respect to 
Fig. |] are noteworthy: Isovector pairing dominates in 
both T — and T — 1 states at low excitation ener- 
gies. With the exception of the states below 1 MeV, 
isovector pairing decreases, approximately exponentially, 
in T = states with energy. On the other hand, the 
isoscalar pairing in T = states is about constant up 
to E — 6 MeV and then decreases more slowly than 
isovector pairing. We conclude that isovector pairing in 
T = states is more concentrated in the low-energy part 
of the spectrum, while isoscalar pairing is also present at 
higher energies. The difference is, as stated earlier and 
shown in detail below, presumably due to spin-orbit split- 
ting, which introduces an additional energy scale into the 
isoscalar pn pairing mode (ej /2 — £5/2 ~ 6 MeV) which is 
only overcome at higher energies. 

The lowest T — 1 states are dominated by isovector pn 
pairing. In these states, like-particle pairing is roughly 
constant, in fact it increases slightly with energy. We con- 
clude that the isovector pairing gap in 46 V is due to pn 
pairing. For states with E* > 4 MeV, when the "pairing 
gap" is overcome, we find no large difference between the 
three isovector pairing correlation strengths. While the 
sum of isovector pairing is still larger than the isoscalar 
pairing strength, the latter is stronger than each of the 
three individual channels. As for T = states, isoscalar 
pairing decreases more slowly than isovector pairing and 
becomes the largest correlation for E* > 15 MeV (how- 
ever, at these high excitation energies our model space 
is quantitatively not adequate anymore and should be 
extended to include the two neighboring major shells). 

What is the origin of the pairing correlations? It is 
wellknown that a "realistic" shell model interaction can 
be approximated by a monopole term and a two-body 
piece consisting of pairing and quadrupole-quadrupole 
(QQ) terms (Bes-Sorensen interaction |25[). Thus the 
correlation can reflect both genuine pairing interaction 
or deformation which is related to the QQ interaction. 
To distinguish between these different sources we define 
first the total correlation energy H corr by subtracting the 
monopole contribution that includes both single particle 
energies and average two body matrix elements |2(f | from 
the expectation value of the full Hamiltonian. Next, we 
measure the importance of the two dominating pieces in 
the residual interaction by calculating the expectation 
values of the pairing, i/ pa ir, and quadrupole-quadrupole, 
Hqq, hamiltonians in the eigenstates of the full Hamil- 
tonian. These Hamiltonians arc defined by 



t z 

-GpioA/ioO) 



(9b) 



where Q M is the mass quadrupole operator defined as: 
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The values of the different constants are taken from ta- 
ble II of reference [p7| . After scaling to the corresponding 
Tiuj the values are: G p0 i = 0.37 MeV, G pl0 = 0.57 MeV, 
and x = 0.014 MeV. The expectation values are shown in 
Fig. 5 as a function of temperature. We have additionally 
split them into their isovector and isoscalar components 
defined by taking the isovector and isoscalar pieces of 
the respective Hamiltonians when they are written in the 
normal order. The results support the discussion given 
above. The isovector correlation energy is stronger than 
the isoscalar one. At the same time, while the isoscalar 
correlation energy decreases gradually and slowly, there 
is a noticeable drop in the isovector correlation energy at 
low temperatures which can be traced back to the pair- 
ing. In the ground state (temperature T = 0) most of the 
correlation energy is due to pairing energy, but the ex- 
pectation value of the isovector component of the pairing 
Hamiltonian has dropped to about half of its ground state 
value already at temperature T = 0.25 MeV. By splitting 
Hpair into its various components, we verify that most of 
the isovector pairing energy (about 60%) and its drop 
with temperature is caused by proton-neutron pairing. 

Strikingly the isoscalar pairing correlation energy de- 
creases very slowly. While the results for the isovector 
pairing correlation energy confirms the expectation and 
shows that this degree of freedom is mainly concentrated 
in states at low excitation energies, our calculation sug- 
gests that a well pronounced "pairing gap" does not ex- 
ist in isoscalar pairing. This conclusion is based on the 
comparison of the isovector and isoscalar pair correlation 
energies in Fig. 5; there is no peak at low energies in the 
T = channel. (The apparent minimum of the T = 
level density around 2 MeV in Fig. 3 appears to be a 
statistical fluctuation.) 

On the other hand, the isoscalar component of Hqq 
increases slightly at small temperatures, counterbalanc- 
ing the slight decrease in the isoscalar piece of H pa i r . The 
isovector component of Hqq is roughly constant in the 
temperature regime T < 0.25 MeV, which is dominated 
by the states with large isovector pairing contributions. 
We also notice that the isoscalar quadrupole-quadrupole 
correlations decrease more slowly than the isovector ones. 

In Figure |^ we compare the values of the different cor- 
relations obtained by adding the isoscalar and isovector 
parts. Clearly pairing plus quadrupole accounts for most 
of the total correlation energy as measured by H corr . 
The contributions of the other components of the nuclear 
force (see |27| for a complete characterization) are prac- 
tically independent of the temperature. We note that 
at T w 0.15 MeV, where the peak in the specific heat 
in figure ^ appears, the quadrupole correlations became 
dominant over pairing correlations as the low lying T = 
states, that dominate the thermal average at this tem- 
peratures, have bigger quadrupole energies than pairing 
energies. 



Finally, we like to report on SMMC calculations of 
the even Fe-nuclei 50_58 Fe. We have chosen this set of 
isotopes as it includes an N = Z nucleus ( 52 Fe), a closed- 
shell nucleus ( 54 Fe), and nuclei with valence protons and 
neutrons in different subshells ( 56 ' 58 Fe). Our calculation 
parallels the SMMC study in Ref. (|] where the reader 
might find details about the approach. Parts of the re- 
sults concerning pairing have already been presented in 
Ref. although in different context. 

To distinguish between genuine pair correlations and 
those reflecting the different number of neutrons, we in- 
troduce the excess of pair correlations M exc defined by 
subtracting the mean-field values from the pair correla- 
tion strength defined in Eq. (^j). Following Ref. |nj we 
define the mean-field value by the uncorrelated Fermi gas 
value using the factorization 

( a I a /3 a 7 a <5) = nffn a (8^S a s - bps^u-i) (10) 

where n a are the occupation numbers. Fig. [7] shows 
the excess of isovector and isoscalar pair correlations. 
Furthermore we have indicated which orbital coupling 
mainly contribute to these excesses. 

Isovector pn correlations are strongest in the N = Z 
nucleus 52 Fe and decrease rapidly with increasing \ N — Z\. 
This result confirms the general trend already outlined 
in jjj which noticed that with increasing \N — Z\ nu- 
clear ground states show the tendency to split into sep- 
arate proton and neutron condensates. The decreasing 
pn correlations also allow an increase in pp correlation 
with growing |7V — Z\, as can be also seen in Fig. 7. 
The excess of nn correlations, however, does not follow 
the simple SO (5) picture as shell effects, which are not 
present in that model, play an essential role. One clearly 
observes that Af exc decreases strongly towards the magic 
neutron number N — 28 reflecting the closure of the fa 12 
shell in 54 Fe. With two neutrons outside of fa 12 neutron 
correlations increase again ( 56 Fe), but are somewhat re- 
duced again in 58 Fe due to the partly closure of the p^/2 
subshell. For this nucleus most of the nn correlations are 
due to pairing in the fa/2 and Pxn orbitals. 

A closer inspection of the pp correlations indicates that 
the /g/2 orbitals contribute more than the p^/2 orbitals, 
although the latter are energetically favored. This is al- 
ready a signal for the presence of the isoscalar pn pair- 
ing which introduces 77/2 — fa/2 correlations which are 
roughly the same for all investigated Fe- isotopes. As long 
as the fa / 2 neutron shell is not filled, we also find isoscalar 
pn correlations among the fa/2 orbitals. For the heavier 
isotopes these correlations are actually smaller than the 
mean-field values indicating that other correlations be- 
come more important. 

As suggested above, isoscalar pn correlations are ex- 
pected to be sensitive to the amount of spin-orbit split- 
ting. To quantify this statement we have repeated our 
SMMC studies of 50 Fe and 56 Fe by artificially reducing 
the spin-orbit splitting in the single-particle energies by 
a factor of 2; thus we use £7/2 = 0.,e3/2 = 2.,e 5 / 2 = 
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3-7 £i/2 = 3. As expected, the reduction of the spin-orbit 
splitting generally increases the isoscalar correlations at 
the expense of the isovector correlations. In particular, 
the isoscalar correlations due to the coupling of the spin- 
orbit partners fj/^ — /s/2 are strongly increased (by about 
a factor of 2.5 in these nuclei). However, a deviation 
from this general trend is also noteworthy. The increase 
of isoscalar correlations leads also to an increase of nn 
correlations in 56 Fe due to increase of correlations in the 
fb/2 orbital. 

Obviously coupling between /V/2/5/2 is very strong. 
Will it dominate if the spin-orbit splitting were removed? 
Unfortunately such a calculation cannot be performed 
with the SMMC method due to the breakdown of the 
g-extrapolation [^0| . For that reason we have gone back 
to 46 V and have performed one shell model diagonaliza- 
tion study with the KB3 interaction and one in which 
the monopole terms have been removed from the inter- 
action. In both cases we then calculated the expecta- 
tion value of the isovector and isoscalar pairing Hamil- 
tonians in the lowest T = state, which is the ground 
state in this academic case. For the eigenstate of the 
realistic Hamiltonian we find (Hj^) = 2.18 MeV and 
(Hpcri®) — 1.41 MeV, while the magnitude is inverted if 
the monopole terms are taken out: (H^*) = 2.09 MeV 
and (Hjcrir) = 4.62 MeV. It is thus obvious that isovector 
correlations in odd-odd N — Z nuclei win over isoscalar 
correlations due to the presence of the spin-orbit split- 
ting. 



III. CONCLUSIONS 

Modern diagonalization codes make it possible to 
study nuclear spectra state-by-state even for nuclei rather 
far away from closed shells, where correlations play a de- 
cisive role. Here, we use the N = Z odd-odd nucleus 46 V 
as a case study for the investigation of pairing and other 
correlations. We are able to describe not only the ground 
state, but the development of various correlations with 
the excitation energy and, respectively, temperature. 

Our main emphasis is on the competition between the 
two basic modes (like-particle {pp and nn) and proton- 
neutron (pn)) of the isovector pairing on one hand, and 
the isoscalar pairing on the other hand. These features 
are particularly important in 46 V, where the ground state 
spin and isospin J, T = 0, 1 demonstrates the importance 
of the np isovector pairing. We show that this mode of 
pairing peaks sharply at low energies. At energies or tem- 
peratures above the isovector pairing gap (E > 3 MeV 
or T > 0.5 MeV) the three isovector pairing strengths 
become essentially equal to each other, and gradually 
decrease. At the same time, the isoscalar pairing energy, 
which is smaller than the isovector one inside the gap, 
becomes at higher energies comparable with the total 
isovector pairing, or even bigger. 

We have shown (Fig. 2) that the large statistical 



weight of the low-lying excited isospin T = states causes 
a characteristic maximum in the specific heat Cy . This 
feature represents an analog of the phase transition from 
the T — 1 ground state to the regime where the partition 
function and internal energy are dominated by the T = 
states. 

Finally, we study the contributions to the correlation 
energy at different temperatures in 46 V. We show, first of 
all, that as expected pairing and quadrupole correlations 
account for most of the total correlation energy every- 
where. All other correlations contribute about 2 MeV, 
and this amount is essentially independent of tempera- 
ture. Pairing dominates near the ground state, but soon 
(for temperatures above 250 keV) pairing and quadrupole 
correlation energies become comparable. 

In an additional study based on the SMMC we show 
how pair correlations in Fe isotopes depend on | N— Z\. In 
particular, we were able to trace the contribution of indi- 
vidual subshells to the pairing strength. That calculation 
also confirmed that the isovector pn pairing strength de- 
creases fast with increasing \N — Z\. On the other hand, 
the isoscalar pn pairing in this chain of isotopes is almost 
independent on \N — Z\. 

In the SU(4) limit, i.e. for degenerate single particle 
levels, identical coupling constant of the isoscalar and 
isovector pairing force and absence of other components 
in the residual interaction, the lowest T = and T = 1 
states are degenerate. However, these symmetries are 
lifted in realistic cases and, for the present nucleus 46 V, 
the ground state has T — 1. At first we note that the 
coupling constant of the isoscalar and isovector pairing 
force, determined in Ref. ^?J, are G p oi = 0.37 MeV, and 
Gpio = 0.57 MeV, i.e. the isoscalar pairing is stronger. 
Why does the isovector pairing nevertheless dominate the 
ground state? We argued that the effect of the isoscalar 
pairing is weakened by the spin-orbit splitting which in- 
troduces additional energy scale into the problem. The 
isoscalar pairing can act efficiently only at energies where 
the spin-orbit splitting plays only a minor role. To prove 
that point we have artificially reduced the spin orbit split- 
ting in Fe isotopes, and indeed observed an increase of 
the isoscalar pairing. But through higher order effects 
the isovector pairing was affected as well. Furthermore, 
we have shown that in the lowest T = state in 46 V 
the isoscalar correlations would prevail if the spin-orbit 
force were switched off. Thus, this somewhat academic 
exercise has confirmed to us that without the strong spin- 
orbit force the nuclear ground states would look very dif- 
ferently from the ones we are familiar with. 
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FIG. 1. Calculated and experimental spectra of V. Only 
the isospin of the ground state is known experimentally. The 
energies of the excited T=l states in 46 V can be estimated 
from the isobaric analog states in 46 Ti: 2+ (0.889 MeV), 4 + 
(2.010 MeV), 0+ (2.611 MeV), 2+ (2.962 MeV). 
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FIG. 2. Internal energy U (upper panel) and specific heat 
Cv (lower panel) of 46 V as a function of temperature T. 
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3. Level density of 46 V as function of the excita- 



tion energy. The levels with isospin T=Q and 1 are shown 
separately. The Fermi gas level density parameters were 
a = 3.5 MeV -1 and A = -2.0 MeV. 
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FIG. 4. Average isovector T=l,S = 0,L = and 
isoscalar T = 0, S = 1,L = pair correlation strengths in 
46 V as function of the excitation energy (binned for clarity in 
bins 0.5 MeV wide). 




FIG. 5. Expectation values of the total correlation 
energy, isoscalar and isovector parts of pairing, and 
quadrupole-quadrupole energies as a function of the temper- 
ature for 46 V. 
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FIG. 6. Expectation values 

of the pairing, quadrupole-quadrupole and total correlation 
energies as a function of the temperature for 46 V. 
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FIG. 7. Isovector and isoscalar pairing correlation in Fe iso- 
topes. The excess over the mean field value is plotted. The 
unfilled symbols refer to the SMMC studies of 50 Fe and 56 Fe, 
in which we have artificially halfed the spin-orbit splitting in 
the single particle energies (see text). 
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